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in 

' Introduction 

-5, 

The path integral method was proposed by R. Feynman as a new approach [1| to the solution 
of quantum mechanical problems. Nowadays it became one of the most powerful methods in the- 
oretical physics. Many applications (see [2H5|) °f this method are devoted to diverse important 
problems. In articles [g, [7( describing physical analogies, the path integral method was applied to 
a series of financial tasks. One of these tasks is the problem of evaluating the option price, initi- 
C-T) ' ated in the seminal work [8j which is known as the Black-Scholes model (see appendix [X]) . Some 

generalization of this model is the Merton-Garman model, where there was suggested ](| 0, @| a 
suitable dynamics equation of the option price. Just like the Schrodinger equation, the Merton- 
Garman equation is of evolution type. Hence, the path integral method is well fit for presentingthe 
corresponding solution in a closed form and reduces the problem to quadratures. In articles [a, @| 
the path integral representation was constructed for the operator evolution kernel (propagator) of 
a Merton-Garman Hamiltonian. 

Meanwhile, while constructing the evolution operator kernel in 0, 0|, there were encountered 
a few shortcomings. In particular, the operator structure of the Merton-Garman Hamiltonian was 
not taken into account which holds the terms with mutually non-commuting summands by analogy 
with momentum and coordinate operators of quantum mechanics. In the case under regard, the 
evolution operator kernel for the Merton-Garman equation is not determined within the Feyn- 
man path integral (on the phase space) from the classical Hamiltonian function entailing from 
the Hamiltonian operator by replacing the momentum and coordinate operators by their classical 
expressions p y — > p y (p y = —id/dy). Correspondingly, the same holds for the path integral in the 
configuration space, and thus the classical Lagrange function does not make it possible to correctly 
determine the evolution operator kernel. The situation mentioned above is well known 

@, SUB-ELI 

in quantum mechanics and concerns the operator ordering problem. It is known (Io| that in the 
case of terms with non-commuting factors within the Hamiltonian the corresponding action in the 
path integral still contains additional components. Subject to the Trotter formula approach (which 
was used in @, 0|), the above mentioned components were chosen from the correct form of the 
resulting Hamiltonian evolution equation based on the obtained evolution operator kernel. As it 
follows from (see also [Hj]), since these additional components were not taken into account the 
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evolution operator kernel does not correspond to the initial one but rather to some other operator 
symmetrized over the momentum operators with the Hamiltonian operator. Thus, the evolution 
operator kernel constructed in does not correspond to the initial Merton-Garman Hamilto- 

nian but to some other Hamiltonian. The determination of these additional components within the 
Trotter approach is a rather cumbersome procedure [10]. In our work, we use an alternative ap- 
proach for construction of the path integral evolution o per ator kernel representation based on the 
Gaussian path integral within the T- product approach [111 Il3| (see appendix [Bj . It is well known 



that this approach is effective for Hamiltonians with quadratic momentum dependence, that is the 
case for the Merton-Garman Hamiltonian. We constructed the corresponding path integral evo- 
lution kernel representation for the Merton-Garman equation having interpreted the integration 
variable as the velocity pi. fll|. Moreover, we determined the additional components within the 
natural Gauss integration procedure assuming the proper symmetrization. Based on this evolution 
operator kernel we presented an analog of the Black-Scholes formula for a European call option 
price. 

1. The basic Merton-Garman model setting 

Within the Black-Scholes model [H, Q two main financial instruments are considered: bank 
account and stocks. A bank account change H(t) > - on risk-free instruments and governed by 
the equation: 

dll(f) = rn(t)di, (1) 

where the value of r ^ - interest rate. A share price change S(t) ^ - can be determined from 
the stochastic equation: 

dS(t) = <pS(t)dt + aS(t)R(t)dt, (2) 

cj) - is the expected return on the security S(t) and a - stock price volatility (the standard 
deviation). The quantities <f> 1 a are taken to be constant. The second summand (0) contains a 
Gaussian random function R(t) (the white noise), whose first moments are equal: 

(R(t)}=0, (R{t)R(t r )) =6(t-t'). (3) 

The brackets above (. . .) denote the Gauss average of the random variable R(t). Stochastic equa- 
tion ([2]) is called a geometric (as well as economical) Brown motion [3] . As concerns the classical 
Langevin equation for the Brown particle velocity [TBI . [l6| , a random component of equation @ is 
proportional to S(t). The quantity dW(t) — R(t)dt in equation © is also referred to as the stan- 
dard Brown motion (or the Wiener process) which describes the independent increases of quantity 
R(t) with characteristics: 

(W(t)} = 0, {W{t)W{t')) = min(i, t'). 

The basic assumptions of the classical option pricing theory are that the option price C(t, S{t)) 
at time t is a continuous function of time and its underlying asset's price S(t). To realize this 
dependence, it is necessary to impose some relationships. The Black and Scholes idea, concerning 
the determination of a rational option price was based on a call option payment (premium) 
C(i) at a moment of time t (0 ^5 t ^ T), considered by an option seller as a minimal capital 
(investment portfolio), whose investments will permit a contract to get fulfilled. Making use of this 
aspect, Black and Scholes have determined d, gj the following investment portfolio: 

C(t,S{t))=R(t) + dC ^ )) S(t). (4) 

The first summand in ^ means the premium part which is placed in the bank account, the second 
part is assigned for buying shares at market prices S(t). Within the economical sense, namely 
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coefficient dC(t, S(t))/dS(t) > determines the amount of shares bought at the price S(t). For 
option price dynamics equation to be obtained we need to find the change in time of both parts 
of equality (H}, having taken into account equations JT]), ^j. In particular, in the left hand side 
of (UJ) one obtains (for brevity, the ^-argument is omitted): 

dC_dC dCdS 1 d 2 C (dS) 2 

dt ~ dt + dS dt + 2 dS 2 dt ' ( ' 

The last summand in ([5]) needs to take [H,[l4j the relationship (dS*) 2 = a 2 S 2 dt into account. Taking 
now the derivative from the right hand side of Q and equating it to (JSJ), upon some simplifications, 
one obtains 0, @, [3| the Black-Scholes equation: 

dc dc i 2 2 d 2 c 

^ +rS dS + 2 a S dS^^ ra (6) 

Making change of variables S — e x (x € K) we obtain equation (03]) . Its solution (appendix |X| 
permits to obtain the Black-Scholes formula for the rational option price (|55p. 

As it was noticed above, the volatility a 2 within the Black-Scholes model is a constant value. 
Thus, some generalizations of the Black-Scholes model were proposed in which the volatility changes 
with time. In particular, in a rather general form, the volatility dynamics equation was written 
d, [§| as follows: 

dV(t) = [A + fiV(t) + £,V a {t)R 2 (t)] dt, (7) 

where V(t) = a 2 (t), i?2(i) _ Gauss white noise. There are also conditions imposed on the constants 
A, fi , £ which provide the positivity V(t) > 0. Then, the dynamics share price can be written as 
follows: 



dS(t) = <f>S{t)dt + y/V(fjS(t)Rx(t)dt. (8) 

The correlation Gauss white noise Ri(t), i?2(£) moments are determined by relationships ([3]), as 
well as by: 

(Ri(t)R 2 (t')) = pS(t-t') (9) 

with the correlation parameter — 1 p 1 . 

While deriving the dynamics option price equation the Black-Scholes portfolio property (jU) was 
taken into account. It is easy to see that the relationship structure of ((4]) gives rise to the summand 
contraction of ^ dS(t)/dt, and the equation © does not contain the terms proportional to the 
random component ^ R(t). Meanwhile, the price S(t) is determined at time t (similar to the Brown 
particle velocity, whose velocity is known at initial moment t, but its change is a random variable). 
Making use of this property while forming the investment portfolio, the following equation was 
obtained (its more detailed analysis can be found in [y, |7|) for the option price dynamics: 

£ + 4§ + ,A + „V)£ + iwg + + _ rC . m 

Introducing new variables S = c x , V = c y (x, y 6 ffi.) equation (fTD|) can be written down as follows: 

dC(t) 



dt 



= H MG C(t), (11) 



where we denoted the quantity: 



2 8x 2 + \2 J 8x + \ + M 2 J dy 





d 




a — h r 
ox 


a 2 




dxdy 


2 



p^"" 1 / 2 ) _iL_ _ ^e 2 ^"" 1 ) |^ , (12) 

which means exactly the Merton-Garman Hamiltonian. It is easy to see that the first three sum- 
mands in p"2")) correspond to the Black-Scholes Hamiltonian term, while the next ones are owing 
to the stochastic dynamics of the volatility. 
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2. The path integral representation of the Merton-Garman equation 



The continuous solutions to equation (|46|) . which is understood in the standard generalized 
sense, can be represented in the form: 

C(t) =c- tHmg C(T), (13) 

where C(T) is the value of the option price at time t, which is determined by function (I5U)) . The 
evolution operator kernel in (|13p . similar to (|48p . can be determined by the relationship: 

g{r, x-x Q ,y, y ) = e~ THMG S(x - x )S(y ~ y ). (14) 

Since the Hamiltonian coefficients in (I12[) are independent of the variable x, based on the 
experience at the construction of the Black-Scholes Hamiltonian and its evolution operator kernel, 
one can make the Fourier transformation with respect to the x- variable and as a result, one obtains 
the Hamiltonian (Ti~2")) in the form: 

H MG (k)^h (k,y) + h 1 (k,y)^-^^ a - 1 ^-^j , (15) 
where we have denoted: 

ho(k,y) = ^k 2 + i(^-r)k + r, (16) 



hi(k,y) = -ip^-Wk - (\c-y + n - i^Lfftl"- 1 )} . (17) 

Respectively, the Fourier image of the kernel (fl"4l is determined as: 

g(r, k, y, y ) = e" r ^ - y ). (18) 

2.1. The case a — 1 

Consider the simplest case a = 1 where the coefficient at the operator d 2 jdy 2 does not depend 
on the variable y. Having chosen the complete square form in the Hamiltonian (|15p . one obtains: 



Hu G (k) = -\ (e| - h ^f) + h (k,y) + - \h\Ky). (19) 



Here the quantity h(k, y) equals the value h\(k, y) (|16p for a = 1: 

h(k, y) = -ip£,e y/2 k - (Xe~ y +p- £ 2 /2) 
and the summand h'(k,y) means the derivative of the h(k,y) with respect to the variable y: 

h'(k,y) = Xe-y -ifcip^ 72 . 

Namely, a similar summand, as we will show below, was not taken into account in articles [6j, [7j. In 
our case, this appears to have been caused by the corresponding symmetrization performed in (|19l) . 
One can also observe that this summand is important for the Hamiltonian (1191) (in the coordinate 
representation k — > —\d/dx) whose form would take an incomplete form: 

e y d 2 f e y \ d ( e\ d x f _ d d _ v 

^fe^ + ^ 2 )-^- (20) 



2 dx \ dy dy J 2 dy 
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It is seen that that summands with non-commuting factors in (|20p are written in a symmetrized 
form. Namely, as it was already mentioned above, the path integral representation of the evolution 
operator kernel, constructed in [6, 7], corresponds to ([201) in contrast to the initial Hamiltonian (fT2")l 
(for a = 1). 

The evolution operator in (|T5jl can be transformed by using the expression (|5()1) : 



7V(r)/ Duexp -- / v 2 (t)dt 



x Texp 



1 



1 



1 



-hoik, y) - 7^h 2 (k, y) + ^h'(k, y) + ^v{t)h(k, y) 



T 

dt- J £v(t)dt 



d_ 

dy 



(21) 



The latter summand in the T-exponent contains a shift-operator which can be disentangled with 
respect to (|63[) . since the T-exponent becomes a usual one: 



Texp - J &(t)dt-^ = exp - J ^v(t)dt^- 



Texp I - J ^v(t)dt^- J = exp I J £v(t)dt 



d_ 
dy I ' 



Respectively, for the operator kernel (|18p . upon simple transformations which are omitted, one can 
obtain: 



g(T,k,y,y ) =Af(r) Vvexp 



1 



[v(t)-h(k,y)] 2 dt- / h {k,y)dt x 



x exp ( - / h'(k, y)dt \S\y-y - / v(t)dt ) , 



(22) 



where we denoted y = y — v(t')dt' and the constant Af(r) is determined from the condition: 



M{t)\ Vvexp ( - — / v 2 (t)dt\ =1 



The path integral (|22l) is defined in the space of velocities 0, [H[ • The quantity 

L(k) = ^[v-h(k,y)} 2 + h a (k,y) 

means a Lagrangian function (for imaginary time) and the summand —h'(k,y)/2 - the additional 
term. Relationship (|22D determines the Fourier image with respect to a;-variable. For the Fourier 
inverse transform one needs to separate in (|22[) the fc-dependence: 



l{r,k,y,y ) = e 



r Af(r) J Vve s ° exp {-\^d 2 k 2 - irk [a 2 /2 - f]\ 5 iy - y - J v 



(t)dt , 



(23) 
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where we denoted 

T 

1 



So 



2e 



[v(t) + Xe-« + (i - /2] 2 dt + -A / e _fi di, 



1 



a 2 = / e»di, 



(24) 



Performing the inverse Fourier transform we obtain the evolution operator kernel: 



g(r, x-x ,y, y ) = 7V(r) / T>ve 



-.So 



V27 



x exp 



1 (x-xq+t [f-a 2 /2])' 



* | //-//„- J viDdl | . (2.-.I 



Thus, formula (|25|) is a path integral representation of the evolution operator kernel for the Merton- 
Garman Hamiltonian in the case a = 1. As it was mentioned earlier, the comparison of this result 
with the formula from [(| 0| shows that the latter does not contain the second summand in Sq , as 
well as there is a change in the corresponding numerical coefficient in the third summand f. Recall 
that in 0, 0] there is used a Feynman path integral defined in the configurational phase space, 
in contrast to the kernel (|25[) which is represented in the velocity phase space. The connection 
between both representations can be found in appendix [Cl 

Making use of this kernel, we obtain a generalization of the Black-Scholes result (|55l) for the 
option price. For this to be done, we need to multiply the kernel g(r, x — xo , y, ya) by the payment 
function (|53[) and integrate over the variables xq, yo- Since function (|53[) does not depend on yo , 
the integral with respect to yo is easily performed owing to the presence of ^-function. As a result 
we obtain: 



C(t,S(t)) =M{t) Vve bo S(t)e TAr N(d + ) 



"KN(d-] 



(26) 



where we denoted Ar = r — r and quantities d± are determined by (|55[) for d± by means of the 
change r — > f, a — > a. 



2.2. The case of arbitrary a-values 

In this case, in order to obtain the path integral representation of the evolution operator we 
will act the following way. First, we will make use of the path integral (|56[) for linearization of the 
third summand in (1151) : 



-T%G(*i) 



Af(r) / Vw exp 



x Texp 



w 2 (t)dt x 



h (k,y)dt - 
o o 



h 1 {k,y)+£w(t)e v ( a - 1) 




(27) 
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Then, we use the path integral (|61j) for linearization of the summand with the derivation opera- 
tor (|2"Tj) . Preliminarily, let us denote: 



l r 



A(t) = -- h 1 (k,y)+tw(t)^ a - 1 '> 



dy 



l r. 



[A(t),B(i)]_ = - [h[(k, y) + (a- l)^)^" 1 )] = 

The sense of the previous transformations, as we observed before, is in splitting a shift-operator 
which permits to perform the T-operation. As a result, for the T-exponent (|27[) we obtain the 
relationship: 



Texp 



h (k, y)dt - J [/n (A, y) + Zw(t) e y( a -V 
o 

M{t) 2 J VvVu exp I ~J v 2 (t)dt-i J v(t)u(t)dt 



d_ 

dy 



dt 



xTexp - / [ 



([v(t) + 2iu(t)]A(t) + v(t)-^ + A'(t) + h (k,y)^dt} . (28) 



Disentangling the exponent containing the shift-operator and substituting it into (|27[) . we write 
down the evolution operator in the form: 



rH MG (k) 



I 



= jV(t) / exp ( -- / t/(t)dt | $(w)exp | - / h (k,y)dt | exp | - / u(*)dt^ 







(29) 



where we denoted: 



$(t>) = Af(r) J Vucxp I -i J v(t)u(t)dt J 7V(r) j Vw exp I / w 2 (t)dl j 



x exp 



v(t) + 2iu(t)]A(t)+A'(t)\ 



dt 



(30) 



Let us note here that (~)-barred quantities differ from those, used before, by the change to y — > y, 
where y = y — J v(t')dt'. 

While calculating $(u) it is necessary to perform integration with respect to variables w(t) 
and u(t). Since the mentioned integrals are Gaussian, they can be performed in closed form. Upon 
simple but a bit cumbersome calculations, which are omitted, we obtain: 



x exp 



M{r)j VvHe-^-V exp [ 



b(t)-fei(fc,y)] 2 

e 2y(a-l) 



dt x 



h (k,y)dt+- / h' x {k,y)dt+-{a-l) / [v(t) - h^k, y)] dt exp - v{t)dt 



0_ 

dy 

(31) 



Respectively, for the Fourier representation of the kernel (fT5|) one obtains the relationship: 



g(T,k,y,y ) 



=M{t) Jvv\{ e HK"-i) exp - y «(t), j/) dt J S I y - y - J 



v(t)dt , (32) 
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where we denoted the function: 

L(k,v,y) = ± [v -j£( k _'V )]2 + h (k,y) - ^(k,y) + \{a - 1) [v - h^y)} . (33) 

Let us mention here that the representation (j3"2")) corresponds to the structure [I, [HI, 03 of the path 
integrals subject to the curvilinear coordinates. As is seen from (|32p the kernel is not determined by 
means of a "classical" Lagrange function which is assigned to the Hamiltonian (fT2"| . It is not hard 
to observe that this Lagrange function (for the imaginary time) is determined by the relationship: 

L (k, », y) = ^ ^ ~2 H V iti) ] + fe o(fc, y), (34) 

where we denoted 

h(k, y) = -ifW^-Wh - (^X c -y +n- ^ c 2y(a-i)^ f (35) 

the quantity ho(k, y) is determined by Note also that h(k, y) is different from h±(k, y) (jTTJ) by 

a factor a in the last summand. It is evident that the difference between L{k, v, y) and Lo(k, v, y) 
determines the above mentioned additional terms. The corresponding terms are not taken into 
account in Q for the constructed evolution operator kernel of the Hamiltonian (fT2")l . Note also 
that at the change L(k,v,y) — > Lo(k,v,y) in the obtained kernel will correspond not to 

operator (TT^)) . but to its analog iJ^G > being in some way symmetrized over operators d/dy. The 
structure -ffj^G can be determined, in particular, making use of the approach devised in 



iigy. 

m 

The next step, having split the dependence on the fc-variable, one can obtain: 



Kr, k, y,y )= e~ T W(r) f Vv JJ e ^ a ^e s ° x 

t 

x exp ^— — r<7 2 fc 2 — irk [cr 2 /2 — r\^j 8 [y — ya — J v 



(t)dt , (36) 



where we denoted: 



T 2 



2^2 y e 2 i'( Q - 1 ) 

o oo 

T 

ct 2 =(l-p 2 )i /" e s di, 
o 

r r r 

>i y e*d*-^i J J [v(t)-h(y)] e^/ 2 -^dt, 



■ r — p 



h(y)=-[\e-y + »-^ r e^ a -V). (37) 

Here h'(y) means the derivative of function with respect to a suitable argument. Performing the 
inverse Fourier transform we obtain the evolution operator kernel in the case of arbitrary a- values. 



g(r, x-x ,y, y ) =tf(r) J Vv J] e^^e 



p—rr 



%/27TTt7 2 



: x 



, 2 * 



x exp 1 (*-«> + ' \T- * W ]6 \y-y -l „(t)d* I . (38) 
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Then, the generalization of the Black-Scholes for the option price will take the form: 

C(t,S(t))=AT(T) f VvWe-^-^e 30 \s(t)e TAr N(d+) ~ e- Tr KN(d^)] . (39) 

3. Some comments concerning path integral calculations 

Works 0, 0j U3| present the formulas for the kernel calculation (|65p in the case a = 1, A = 0. 
Meanwhile, the sub-integrand in (p5|) was represented as a three-tuple series: 

1 / 1 (x-xo + t [r~a 2 /2]) 2 \ 



=^= ex P 



£ onmk I 1 - je^dt\ I ± J e^/ 2 di ) (e^ 2 )\ (40) 



n,m.k 



Then, for a general series term: 

oo y / t \ " / r \ m 

dy / I?'ye s « 1 J e»<*>df I ± / e^)/ 2 dt (e^°>/ 2 ) * 

-oo \ / V / 

there was obtained its closed form. Nonetheless one can be easily convinced that such an approach 
to calculations is wrong, since the used series expansion (|4T)]) does not exist. 

Thus, it is reasonable to mention some approximate methods of calculating the path integrals, 
in particular, for the case a = 1. Consider a formula for the option price (|2l)]) . We will represent 
the corresponding path integrals in ((26)) as those with respect to the Gauss measure: 

C(t, S(t))= ( d»(v) [S(t)F+ (v) - e~ Tr K F_ (v)] , (41) 



where d/i(v) — N'(r)'Dve SG - means a probabilistic Gauss integration measure. Here, the quantity 
Sq is given by the expression: 



r 

= [v(t)+n-e/2] 2 dt, 



and the next terms F±(v) hold the remaining factors of formula (f2l?|) . It is easy to see that the 
values F±(v) are bound by v £ R. Therefore, the integrals used in (|41l) exist. For an approxi- 
mate calculation of path integrals (|41l) one can use numerical quadrature formulas, presented, in 
particular, in [l8| . 

The use of approximation approaches is also of some interest. In solving the problems of quan- 
tum mechanics and statistical physics, the mean-path-approximation 0, lig | was effectively used. 
From this point of view let us rewrite formula (|26p in the form: 



C(t, S(t)) = AT(r) j Vv [S(t)e s + - e~ Tr Ke s -] . 

The mean-path approximations could be naturally considered as zero-approximations. In this case 
within the quantities S± we will make the change: 

r r 

y(t) -> y m = - / y(t)dt = y-~ { tv(t)dt. 
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Upon the above mentioned change, the quantities S± will depend on the integration variables via 
the expressions: 



v(t)dt, 



tv(t)dt 



and the path integrals in (|26|) will then reduce to the usual quadratures. As a result, the price 
option formula will take the form: 



C(t,S(t)) 



/2)/2 



dwid«2 

( 

2tt 



S'(i)e TAr N{d* + ) -e~ Tr KN(d*_) , (42) 



where the corresponding quantities are defined by the expressions: 



Sn* = 



±- (2vf + 6v 2 2 - 6v lV2 - TV! [fi C 2 /2]) - ^e-y + ^(n £ 2 /2)- 



A 2 r 



-2(y— u 2 ) 



A 



At 



dl = 



2^ 2 e 2 

In +r(r* ±ct* 2 /2 



+ - ? e- J/ (l-e^) + ^-e 



cr* = (l-p^e^ 2 , 
r* = r - ^p 2 e y ~ V2 + p 
r =r — r. 



Our next step is to take the series corrections into account, if we expand the quantities S± into 
series with respect to suitable deviations Ay(t): 

Jj(t) = Vrn + Ay(t) . 

Certainly, the approaches mentioned above need a special investigation. 



4. Conclusions 

Using the Gaussian path integral, a representation of the evolution operator kernel for the 
Merton-Garman model was obtained. Also, some shortcomings that occurred in the cited sources 0, 
0| were explained and corrected. A generalization of the well-known Black-Scholes formula for 
the rational price of European call option was also presented. Approximate approaches to the 
calculation of path integrals were also described. 

Note also that the Black-Scholes formula (|55|) is widely used in practice in analyzing the op- 
tions [§]. It is known that the analyzed model gives rise to quite good results in the case of well 
developed markets. However, shortcomings occur when a stock market is in an uncertainty state 
and the price dynamics swings [20|]. At the same time, the Merton-Garman model should be consid- 
ered as more realistic generalization of the Black-Scholes model taking into account the arbitrary 
volatility nature (formulas (0 , © ) ■ That is why the analysis of the obtained ((26)) and ([39)) formulas 
presents a considerable interest which will be the subject of our next investigations. 
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A. The Black-Scholes model 

A.1. The Black-Scholes equation 

The Black-Scholes equation 0, [ljj describes the C(t) option price dynamics: 

dC(t) 



H BS C(t), (43) 

where Hbs denotes the following Black-Scholes Hamiltonian: 

a 2 d 2 fa 2 \ 8 

H * s = -Td^ + [Y- r )d- x +r - (44) 

Equation (l4"3")l visually resembles a Schrodinger equation, written with respect to imaginary time, 
or a diffusion equation. Quantities r and a are considered to be constant. The general solution of 
equation (I43[) is represented in the form: 

C(t) =e tffBS C(0), (45) 

where C(t) is the value in time t, C(0) is an initial condition. When using equation (l4"51 in order 
to evaluate the option prices it is convenient to study the dynamics backward in time. For a given 
value of C(T) at the final moment of time T we need to evaluate C(t) at t < T. It follows from 
the solution of equation (j4"5"]) that 

C(t) = e- (T -* )ffBS C(T) = c- tHbs C(T), (46) 

where we denoted r = T — t. Making use of the kernel g(r, x — Xq), for expression (|46[) one obtains 
that: 

oo 

C(t,x)= J g(T 1 x-x )C{T 1 x )dxo. (47) 

— oo 

Based on (14TH) and (|4T)> one finds the relationship: 

g(r,x-x )=e- THBS S(x-x ). (48) 

Since Hamiltonian (|4"4")l coefficients are constants, it is convenient to proceed to Fourier transform 
by the formula: 

oo 

/(t )!C ) = -L / f(t,k)e ikx dk. 



2tt 

— oo 

In equation (|43l) it is convenient to proceed to Fourier transform, receiving the following result: 

dC(t,k) 



Of 



H BS (k)C(t,k), (49) 



where C(t, k) denotes Fourier transform of proper function and Fourier transform of Hamilto- 
nian (|4"4"1) : 



a 2 



H B s(k) = yk 2 +i( y-rjfc + r. (50) 
An equation (I49p for a Fourier transform will be of the form: 

C(t, k) = e- TilBS ^C(T, k) . (51) 
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Respectively for the kernel Fourier image g(r, x — x$) one has: 

g( T ,k) = e - T " Bs{k) . 

Similarly, for the kernel of the evolution operator (|49p . applying inverse Fourier transform, we 
receive the following expression: 

oo 

g(r,x-x ) = ^ J 5(r,fc)e ife (— <*dk = 



OO 

J exp ^-i T( 7 2 fc 2 - irk [cr 2 /2 - r] ^ e^-^dfc 



m i fe-g + Mr-^])' \ (52) 



\f2 r KT<7 2 \ 2 

A.2. A Black-Scholes formula 

Let us recall a classical example of evaluation for the rational price of the European call op- 
tion 0, [3]. A boundary condition (a payoff function) is defined in the form: 

C(T, S) = (S-K)+ = {%- K, S J (53) 

where if is strike price, S 1 is a market price at a moment T (time to maturity). The payoff func- 
tion (|53p has a sufficiently obvious meaning - option transaction realizes if the stock market price 
is greater than the strike price. Option price at time t < T is evaluated making use of the formula: 



C(t,S(t))= I g(ln[S(t)]-x )(e x ° - K) + dx = 

j- T \ r _ 0-2 /2l) 2 " 
exp -- — 5-^ (e x ° - if )dx , (54) 



where 5(£) denotes stock price at time t. 

Having completed the necessary computations we receive, for option price the well known 
Black-Scholes formula: 

C(t, S(t)) = S(t)N{d+) - Ke- rT N(d_). (55) 

We have here denoted: 

N{x) _ _L / e -V, dz , d± . HS(t)m + r (r±,y2) 

27T J (7y/T 

— oo 

B. Gaussian path integrals 

The Gaussian path integral is used for the degree linearization in an T-exponential and looks 
as follows: 

(T T \ 

J^A 2 (t)dt + J B{t)dt 
o o / 

= M(t)J Vvexpi-^J v 2 (t)dt J Texp I - J v(t)A(t)dt + J B(t)dt\ . (56) 
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In this work A(t) and B(t) - differential operators (in general not commutative), "the constant" 
7V(r) is determined from the condition: 

Af(r)[vvexp I ~J v 2 (t)dt\ = 1, Vv = dv(t) . 
\ o / 

The Gaussian path integral (|56[) was used in in constructing the path integral Green function 
representation for some quantum mechanical Hamiltonians. 

It is evident that the field of application of the formula (1551) is wide enough. The path in- 
tegral (|55|) is interpreted as an operator identity. In particular, one can "disentangle" the factor 
Texp (f Q B(t)di) in (|56p in accordance with (|63|) . and obtain the equivalent equality: 

Texp |i j A 2 (t)dt^ =N{t)J Vvcxp v 2 (t)dt^j Texp ^- j v(t)A(t)dt \ , (57) 



where 



A(t) = 



Texp I y B(t')dt' J A(t)Texp ( / fl(t')dt / 



The integration in (|57p is realized with respect to the Gaussian probability measure, since: 

N{t)J Vvexp f~§/ w2 W d< j Tcx P w(*)^(*) d *j = J d/i(v)Texp ^- ^ v(t)A(t)dt \ , 

where d/*(«) =7V(r)X>wexp (-| J T v 2 (i)di). 

To check the identity (|57[) . we expand the left and right hand sides in degree series. In particular, 
T-exponent on the left hand side (|57l) is represented as the degree series 

Texp 1 1 J A 2 (t)dt\ =E^T(iy ^ 2 (t)di ) . (58) 



n=0 



Expanding the right hand side (|ST)) Texp ^— J Q T w(i)A(i)diJ in series and using the Gauss measure 
properties, one can derive expression (|58l) . Really, for the left hand side of equality (1571) one has: 



(T \ CXD T T 

- J v(t)A(t)dt \ =J2 fdtt... J dt„^j 
/ ™=°0 



,-MnT .1!/,.] !(/„.) . IT,!)) 



where /Lt n is the n-th Gauss measure correlation moment. As is well known, all uneven moments 
vanish, and the even ones are expressed by means of the second correlation moment: 



fx 2 = J d/i(u)«(*i)u(t 2 ) = S(h - t 2 ). 
The moment of 2 rt-th order is expressed by means of the moment fj, 2 as the relationship: 

2n = f d/i(«)«(ti) . . . v(t 2n ) =J2^2---^2- (6°) 
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A summand in (|60[) is taken over all possible different smashes of 2n elements by two-element 
subsets (see, for instance, [2l|). The amount of summands in ([50)1 equals 

(2n)! 



2 n n\ 



Making use of symmetry properties (the permutations of operators under the T-product sign) 
it is not hard to demonstrate that all these summands give rise to the same result, thereby the 
formula reduces to ([551) . 

Based on ([55| and the operator identity: 

2A(t)B(t) = \{A{t) + B{t) f - l -{A{t) - B{t)f + [A(t), B(t)]_ , 
let us linearize the product of operators A(t)B(t): 

Texp (j [2A(t)B(t)+D(t)]dt\ = Af(r) 2 J VvVu exp j v 2 (t)dt-^j u 2 {t)dt 

x Texp ^- J {[v(t)+iu{t)]A(t) + [v(t) - iu(t)]B(t) - [A(t),B(i)]_ - £>(*)}dtj . (61) 

Here we have used the following notation of the operator commutator [A(t), B(t)]- = A(t)B(t) — 
B(t)A(t). Using the shift v(t) —} v(t) + iu(t) the right hand side of equation ([6"Tjl can be written 
in the form: 

N{t) 2 J VvVu exp I -i J v 2 (t)dt - i J v(t)u(t)dt J x 
V o o / 

x Texp j - j {[v(t)+2iu(t)]A(t)+v(t)B(t) - [A(t),B(t)]- - D(t)\dt J . (62) 



In the case of differential operators, the relationship (|61[) makes it possible to decrease its differen- 
tial order. Note also that the path integral representation (loTj) is equivalent to the Stratonovich- 
Hubbard identity [22| for the density matrix in quantum statistics problems. 

Jointly with the introduced formulas, the following expression is used to extract T-exponent 
operators [13]: 

Texp(j [A(t) + B(t)]dt \ = Texp(j A(t)dtj Texp(j B(t)dtj , (63) 

where 

/ * 

B(t) 



TexpyJ A(t')dt'j B(t)Tcxj) yj A(t')dt 



t 



C. A connection between path integral approaches 

Consider the case for parameter values a = 1, A = 0. Based on (PJ51 one can write down the 
kernel expression as: 



I , <? e- rr / 1 (x-x + t \f-a 2 /2]) \ , s 

g(T,x-x ,y,y )= / V'ye s °^== = exp --^ °- L . 64 

yo v ' 
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The measure of integration in (l64|) satisfies the condition: 



J V'y exp 



1 

2^2 



-> $(y-yo), 



at the boundary values y(0) = yo , y(r) = y. Moreover, in the above expressions (I24j) one should 
make the following substitutions: 



y->y(t), v(t)->y(t). 



Thereby, we obtain: 



oo y 



g{r,x-x ,y)= I dy V ye 



V 27TTCT 2 



exp 



1 (x-a;o + r [f-CT 2 /2])' 



where 



o 

r 

(1- P 2 )i / e»Wdt, 



a 2 



(65) 



(66) 
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KoHTMHya/ibHe npeflCTaB/iem-m flflpa onepaTopa eBO.moi4iY b Mop,en\ 
MepTOHa-KapMaHa 

Jl.O. G/iawi/ieBCbKi/itP, B.C. ^HiLueBCbKi/iiP 

Kacpeflpa TeopeTWHHoT 4>i3HKH, JlbBiBCbKi/ii/i Hai4ioHa/ibHHH ymBepcwTeT iwieHi iBaHa Oparno, 
Byn. flparowiaHOBa, 12, 79005 JlbBiB, YKpaTHa 

Kacpeflpa eKOHOMinHoT Ki6epHeTMKW Ta iHHOBaTH kh , flporo6ni4bKHti flepx<aBHHw neflaroriHHUw ymBepci/iTeT 
iMeHi IBaHa <t>paHKa, Byn. /led yitpaTHKM, 82100 flporo6i/in, YKpaiHa 

B MeTOfli kohtm Hya/ibHoro iHTerpyBaHHn no6yflOBaHO flflpo onepaTopa eBoruouJi p/\n raMi/ibTOHiaHy MepTOHa- 
KapwiaHa. Ha ocHOBi flflpa OTpwMaHa cpopMyjia p/in h i/i ormioHy, U40 y3ara/ibHtoe BiflOMy cpopwiyjiy B/ieKa- 
LUoy/ica. BKa3aHO Ta ko>k Ha mo>kjim Bi cnoco6w Ha6/ii/i>KeHoro oGni/icneHHH kohtm Hya/ibHi/ix iHTerpa/iiB. 

K/iKJHOBi c/ioBa: KOHTWHyanbHe iHTerpyBaHHU, nflpo onepaTopa eBon>ou,iT, ujHoyTBopeHHH onu,\oHiB, 
cpopMy/ia BneKa-llioynca, Mop,enb MeproHa-KapMaHa. 
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